Response functions of cold neutron matter: density fluctuations 
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We compute the finite temperature density response function of nonrelativistic cold fermions with 
an isotropic condensate. The pair-breaking contribution to the response function evaluated in the 
limit of small three-momentum transfers q within an effective theory which exploits series expansion 
in powers of small q/pF, where pf is the Fermi momentum. The leading order 0(q 2 ) contribution is 
universal and depends only on two fundamental scales, the Fermi energy and the pairing gap. The 
particle-hole Landau Fermi-liquid interaction contributes first at the next-to-leading order 0(q 4 ). 
The scattering contribution to the polarization tensor is nonperturbative (in the above sense) and 
is evaluated numerically. The spectral functions of density fluctuations are constructed and the 
relevance of the q 2 scaling for the pair breaking neutrino emission from neutron stars is discussed. 



I. INTRODUCTION 

The long- wavelength, low-energy dynamics of 
fermionic systems is determined by their response 
functions to soft perturbations, which are characterized 
by length scales that are large compared to the inverse 
Fermi wave vector and energies that are small compared 
to the Fermi energy. At zero temperature the response 
functions of pair-correlated nuclear matter have been 
studied long ago by Larkin and Migdal [l[ . Recently, the 
response functions of pair-correlated nuclear systems at 
nonzero temperature received attention in the context of 
neutral current neutrino emission via pair breaking and 
formation in compact stars and neutrino scattering 
in supernovae [f|. The evaluation of the response func- 
tions involves typically a resummation of infinite number 
of finite temperature ring diagrams. In the unpaired 
limit of normal Fermi liquid, this resummation scheme 
reduced to the familiar random-phase-approximation 
(RPA); for recent applications in nucleonic and neutron 
star matter, see Refs. 

In attractive, cold, fermionic systems the gap A in the 
quasiparticle spectrum is small compared to the Fermi 
energy and the hierarchy of energy scales depends on 
the magnitude of the perturbation, which can take arbi- 
trary values with respect to the pairing gap A. In this 
work we focus on density perturbations and show that 
the two distinct contributions to the response function 
through the scattering and pair breaking processes are 
effective below and above the energy threshold 2A, i.e., 
the energy needed to break a pair. The pair breaking 
processes are of special importance for applications in 
compact stars; to evaluate them, we propose a new sys- 
tematic low-transferred-momentum expansion of the re- 
sponse function which builds on the previous work on po- 
larization tensors of cold superfluid fermionic systems [2j ■ 
Specifically, we show that the pair breaking contribution 
to the polarization tensor posses a well-defined expansion 
with respect to the ratio of the momentum of the external 
current to the Fermi momentum of the fermions. We also 
adopt more general ansatz for the driving terms in the 
integral equations of Ref. [2f by lifting the degeneracy 
among the particle-particle and particle-hole channels. 



We work in the nonrelativistic limit, i.e., the ratio of the 
Fermi velocity to the speed of light is small {vf/c <C 1). 

The density response functions can be utilized to deter- 
mine the spectrum of the collective modes, the stability 
of the system toward clustering, the rates of electromag- 
netic and weak radiation processes, and so on. In partic- 
ular, the rates of neutrino reactions in stellar interiors can 
be expressed through the response function of underly- 
ing matter to vector and axial- vector weak currents [lll - 
1 141 ] . In nonrelativistic limit, the vector and axial- vector 
responses are mapped onto the responses to the den- 
sity and spin-density perturbations, respectively. The 
response functions in the superfluid neutron matter were 
recently computed and the neutrino emission rates were 
determined in Refs. Phenomenologically, these are 

important in modeling the cooling of intermediate age 
neutron stars and the superburst in accreting neutron 
stars [H-Oll- 

It is now well established that at zeroth order (q = 0) 
the pair breaking density response function vanishes, as 
required by the /-sum rule for the polarization tensor, 
which is a direct consequence of the baryon number con- 
servation. Some authors found analytically the leading- 
order contribution to the polarization which arises at or- 
der q 4 [3|, [l8[ . Since there is no general argument that 
requires the coefficient of the leading-order q 2 term to be 
zero, the answer may depend on the approximations in- 
volved in the theory. In section UlI Al we study in detail a 
new small q expansion of the polarization tensor obtained 
in Ref. [2| and find that the coefficient of the q 2 term in 
the series expansion of the density response function is 
indeed nonzero. Furthermore, it turns out to be univer- 
sal, i.e., it depends only on two fundamental scales of 
the problem, the Fermi energy and pairing gap, and is 
independent of the strength of particle-hole interaction. 

This paper is organized as follows. In Sec. [TT| we de- 
rive the vertex functions and polarization tensor in a 
more general setting that in Ref. [2[ by using different 
particle-particle and particle-hole interactions and clar- 
ify the approximations that arise in the weak coupling 
Bardeen-Cooper-Schrieffer (BCS) limit. A small momen- 
tum transfer expansion is applied to the pair breaking po- 
larization tensor in Sec. IIII Al We also show in Sec. MI Bl 
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the results of exact numerical evaluation of the scattering 
part of the polarization tensor. Further in Sec. IIII CI we 
verify that the unpaired and uncorrelated limits are re- 
covered from the scattering part. Our conclusions are 
collected in Sec. |IV] Details of calculations are pre- 
sented in Appendices lAl and iBl We use the natural units 
h = c = 1 throughout and assume that the Boltzmann 
constant Hb = 1. 



II. DENSITY RESPONSE FUNCTION 

In this section we derive the general form of the 
density-response function of neutron matter at nonzero 
temperature. At densities below the saturation density 
neutron matter forms a pair condensate and can be 
described by the weak coupling limit of the BCS theory. 
This affects not only the approximations that are applied 
to the gap equation, but also the approximate relations 
between the loop integrals, as we discuss below. 

The couplings in the particle-particle (pp) and particle- 
hole (ph) channels, v pp and v p h, are assumed zero range; 
often their values are taken to be degenerate v pp = v p h 
and equal to the lowest order Landau parameter /q. Wc 
shall lift this approximation below by assuming v pp 
v p h- The spectrum of paired neutrons is given by 

e p = v /e p 2 + A 2 (p), (1) 

where £ p = p 2 /2m* — pi is the quasiparticle spectrum 
in the normal state, A(p) is the energy gap, p is the 
three-momentum, m* is the effective mass and fi is the 
chemical potential. For contact pairing interaction the 
gap is momentum independent, A(p) = A. 

The softness of the modes implies that their wave vec- 
tor \q\ <C Pf- Accordingly, we write 

c P 2 (, , 2p-g q 2 \ 

2m* \ p l p l J 

and consider the second and third terms in the bracket 
as small compared to unity, since p ~ p F , where pp is 
the Fermi momentum. Thus, we may write — £ p + 
Ho(2yx + y 2 ), where x = (p • <?)/(H|g|) and y = q/p F , 
whereby y -C 1. Here pto = p 2 F /2m* is the chemical 
potential at zero temperature; we shall drop the index 
in the following. Several observations are in order: 

1. If the expansion is carried out with respect to the 
small parameter — £, P + q — £, P , as in Ref. 0, the 
power counting is not manifest. At the leading- 
order the terms which scale linearly in x drop on 
angle integration in symmetrical limits. The only 
nonzero contribution proportional to q 2 is then fur- 
nished by the recoil term. One needs to carry out 
the St; expansion at least to second order to obtain 
all relevant terms that are of order q 2 . 



2. An alternative to the expansion with respect to 
6^ is the expansion with respect to the ratio y — 
\q\/pF- This expansion exploits the softness of the 
modes and applies both for timelike and spacelike 
momentum transfers. It is an alternative to ex- 
pansions in the ratios ui/qv F and qv F /uj, which are 
valid in these regimes, respectively. 

3. Organizing the expansion in powers of vp/c 1 
(nonrelativistic fermions) does not guarantee per se 
the convergence of the series. At any fixed density, 
the Fermi velocity v F is constant and for sufficiently 
large momentum transfers (q > p F /2) the series 
will fail to converge. 

4. Finally, the smallness of the expansion parameter 
is necessary but not sufficient condition for the con- 
vergence of the Taylor series. The validity of the 
expansion should be checked by an exact numerical 
computation of the loop integrals. 

In the following we shall demonstrate that the pair break- 
ing part of the response function can be expanded sys- 
tematically with respect to the y parameter. Such expan- 
sion is thus valid for small three-momentum transfers, 
but arbitrary energy transfers. 

A. Vertex functions 

We start with integral equations for the vertex func- 
tions and derive a (slight) generalization of their coun- 
terparts in Ref. Q that distinguish the particle-particle 
and particle-hole interactions. These we write as sums of 
central and spin-spin interaction terms 

V pp ~ vpp + v pp (a ■*') + ..., (3) 
V ph ~ v ph + v ph (* ■*') + ..., (4) 

where the ellipses stand for the tensor and spin-orbit 
terms that are subdominant at relevant densities in neu- 
tron matter. 

The integral equations defining the scalar vertex, which 
we write in an operator form, are given by [l|, 0] 

l\ - To + VphiGTiG + FtsG + G^F + FTiF), (5) 
f 2 = v pp (GT 2 Gt + Fr 4 Gt + GTiF + FT 3 F), (6) 
T 3 - v pp (G^t 3 G + FTiG + G^T^F + FT 2 F), (7) 
f 4 = T + v ph {G^r 4 G^ +FT 1 F + FT 2 G^ + G^t 3 F). 

(8) 

Here F = —ia y F, <j y is the second component of the 
Pauli matrix, r = 1. When v p „ = v p h, Eqs. ©-([8]) 
become identical to those of Ref. ]2| . Let us now define 
the "elementary loop" as 

-JLx(p)X'(p + q), (9) 
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where X e {G, G\ F, F^} and g is the degeneracy factor, 
which we omit in the intermediate equations and restore 
in the final ones. Direct calculations show that 

n G t F = n F G, ~Kgf = ~Rfgu n GtG t=n GG ; (10) 



these equalities imply that Ti = T4, which is a conse- 
quence of the time-reversal invariance of the system. 



The remaining equations read 



1 - v ph [Tl GG - TL FF } Vp h Tl GF v ph Il FG 
-2v pp U GF [1 - %,n GGt ] v pp U FF 
-2v pp n FG Vp P U FF [1 - v pp n G t G ] 











i-i 












(11) 



There are six distinct loops in Eq. dill) , namely II GG , 
IIff, n GF , n FG , II G t G , n GG t. In the weak coupling 
limit Hgf — — rifGi II G t G ~ n GG t, an approximation 
discussed in detail in Sec. Ill Bl This reduces the number 
of equations from three to two 



where 



1 - v p h A v ph B 



-v vv B 



„C 



r 2 



r 
o 



(12) 



When v p p = v p h these reduce to Eqs. (1G) and (17) of 
Ref . @ • It seen that the interaction in the pp channel is 
absorbed in the gap equation and it is the ph interaction 
that enters the renormalization of the one-loop polariza- 
tion tensor. This result could have been anticipated from 
the limiting form of the RPA polarization tensor of nor- 
mal Fermi liquids (see Sec. IIII Cl) . 



A(q) = n GG (q)-n FF (q), (13) 
B(q) = 2U FG (q), (14) 
C(q) = n GGi (q)+U F F(q)-(vpp)- 1 , (15) 

with q = (ui,q). The solutions for the remaining two 
vertex functions reads 



ri(«) 
r 2 («) 



T C(q) 



C(q) - v ph [A(q)C(q) + B(q)*Y 

r B(q) 

C(q) - v ph [A(q)C(q) + B(qy 



(16) 
(17) 



B. Polarization tensor 



The full polarization tensor is given by Eq. (35) of 
Ref. with the replacement v — > v p h 



U R (q) = 



A(q)C(q) + B{qf 



C(q) - v ph [A(q)C(q) + B(qy 



(18) 
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The "elementary loops" are denned explicitly as 



Tl GG (q) 



n FF ( q ) 



d 3 P 

(2ir) 3 



u 2 p u\ 



2 2 



2 2 

u Pp 



iq + e p - e k iq - e p + e k 

[f(-e P ) f(e k )} 



[/(<*) - f(e k )} 



2 2 



iq-t p - e k iq + e p + e k 



d 3 p 

W) 



2 \ up up 



k 



hp Up 



iq + e p ~ e k iq - e p + e k 

,.2 



[/(<*) - /(<*)] 



iq-e P - etc iq + e p + e k 



d 3 P i i 



1 



[/HpWM^. 
1 



+ 



1 



n GtG(g) = - 



+ 



+ e p + e fc iq - e p - e k 
d 3 p 



iq + e p ~ e k iq - e p + e k 
[f(-e p )-f(e k )]\\ 



[/(e P ) " /(**)] 



(27T) 



2 2 



2 2 



2 2 



iq + e p - e k iq - e p + e k 

[/(-ep) - /(<*)] 



/(e P ) - /(<*)] 



2 2 



iq - e p - e k iq + e p + e k 



(19) 



(20) 



(21) 



(22) 



where k = p + q and the coherence factors are given by 
u 2 p = (1/2) [1 + i p /t p ] and v 2 = 1 - 

The expression for II^ G (g) above applies at arbitrary 
couplings; however Eq. (|14[) presumes weak-coupling ap- 
proximation, because II G j?(<2) = — ILpg(q) holds only 
in this limit. The weak-coupling limit for the func- 
tion HpG{q) is obtained on substituting in square braces 
u 2 — u 2 = v 2 = — 1/2. This statement is equivalent 
to ignoring integrals of the type 



I = 



d 3 p A 2 £ p f(e p ) — f(e k ) 
(27r) 3 2e p e p iq + e p - e k 



(23) 



The integral vanishes because £ p changes sign for mo- 
menta above and below the Fermi momentum, while the 
remainder of the integrand is an even function in the 
vicinity of pf- Thus the loop polarization function re- 
duces to 



IW?) = -77 



1 f d 3 p 



2 J (2*-)3p"> 

1 1 



iq + tp — e k iq - e p + e k 
1 1 



[/(ep) - /(<*)] 



J,p Up 



iq - e p - e k iq + e p + e k 



[/(-e P ) - /(<*)] 



(24) 



The form of the polarization loop ^g^g(q) i s valid for 
arbitrary couplings. The relation n G tG(?) — ^gc(q) i s 



established on noting that, for example, 



W P V . 



p 

u k v p)~^ u k v p anc ^ that the combination in braces vanishes 
in the weak coupling as it leads to an integral of the type 



The contributions to the polarization function due to 
quasiparticle scattering and pair breaking separate. For 
the contributions from scattering the poles are located at 
±(e p — e k ) and the distribution is given by the combina- 
tion f(e p ) — f(e k ). For the pair breaking contributions the 
poles are located at ±(e p + e k ) and the distribution is pro- 
portional 1 — /(e p ) — /(efc). The pair breaking contribu- 
tion vanishes in the limit T — > T~ (T is the temperature, 
T c is the critical temperature of the phase transition). 



Now we write the retarded polarization functions 
H R (q) in terms of the functions A(q), B(q) and C(q) af- 
ter performing analytical continuation (iq — > lu + iS) in 
functions U G G(q), H-FF(q), IIf g(<?) n G t G (<j). After some 
straightforward algebraic transformations we find 
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A(q) = 



B(q) 
C(q) 



d 3 p jl 



(2tt) 3 1 2 



ofa " £ p) 



A 2 



/fa) - /fa) 

{uj + iS) 2 - (e p - e k ) 2 



1 



+ ofa + efc) 



2 

-Aw 



2 1 



f(~e P ) ~ /fa) 1 

(w + Z<5) 2 -( £p + 6 fc ) 2 [> 



d 3 p 1 



/fa) -/fa) 



/(-e p )-/fa) 



(w + i<S) 2 - (e p - e fe ) 2 (w + ^) 2 - fa + e fe ) 2 



d 3 p 



fa - efe) 



x _|pCfc_ A^\ 



/fa) ~ /fa) 
(w + id) 2 - fa - e k ) 2 



fa + e fc ) 1 + 



6>& , A 2 



f(-e p ) - /(e fc ) 
(w + i<5) 2 - fa + e fc ) 



, - [1 - 2/(e P )] 



(25) 
(26) 



(27) 



Note that limg_,.o,w->-oC(g) = 0, since the coupling con- 
stant in the particle-particle channel can be expressed 



as 



{v pp ) 1 



n G tc(9 
d 3 P 



0) 



(2tt) 3 2e, 



U FF (q = 0) 
- 2/fa)] , 



(28) 



which is the gap equation for the contact interaction v pp . 
We do not need to specify the regularization of the gap 
equation, since its divergence is eliminated in the loop 
integrals. For numerical purposes we will adopt gaps 
obtained from finite-range interactions in Ref. [19J. 



III. EVALUATING RESPONSE FUNCTIONS 

Equations (l2"5j) - (j2"T|) separate into scattering and pair 
breaking contributions. We shall see that the first con- 
tributes essentially below the pair breaking threshold 
u < 2A, whereas the second contributes for lj > 2A. 
In the following we discuss in detail the pair breaking 
part and its small momentum expansion. The scatter- 
ing part will be addressed later in this section, where we 
evaluate it numerically. 



A. pair breaking response function: small 
momentum expansion 

The small-?/ expansion of the pair breaking part of the 
polarization tensor is obtained upon writing H R (q) = 
Sn'PnSA A(q) = J2 n -A n y n and similarly for the func- 
tions B and C and truncating the Taylor series at the de- 
sired order in y. The odd powers of y do not contribute 
to the series. 

Up to order y 4 the coefficients for the polarization ten- 



sor are given by 
V = 0, 



V 2 = 
Vi = 



2B0B2+A2C0+A0C2 
Co 

B 2 2 + 2B B i + A i C + A 2 C 2 +A C 4 



Co 



- V 2 



C2 -v 
-7- - v ph V 2 



(29) 
(30) 



(31) 



The coefficients of the expansion are given by (for details 
see Appendix lAf 



An = 2A 2 



Bp, = Aw 



d 3 p 1 

d 3 p 1 



La, 



A 2 = 



Co 
d 3 p 
2A 2 /i£, 



uj 2 f d 3 p 1 



2 J (2tt) 3 e p 



La, 



(32) 
(33) 
(34) 



A 2 /ig P 



(6/x£ p x 2 - e 2 )L 



r 2 2A 2 T 



B 2 
C 2 



(27T) 3 e p 



d 3 p 
(2^)s 



M^p t . 2/i^p 2 1 o t 

— - La H -Li x + 2e p L 2 

£p £p 



(35) 
(36) 

1 

(37) 



where the functions La, L\, and L 2 are defined by Eqs. 
(|A10j) - (|A12| of Appendix El The zeroth order term Va 
vanishes as a consequence of the /-sum rule. The leading 
order nonzero term is given by 



V 2 = A 2 v{p F ){Io + h +h), 



(38) 



where v(pf) = gm*PF /27r 2 is the density of states at the 
Fermi surface, and 

J (27r) J e p uj 2 



» + p( ! f t *"- 1 

j 9 , n s-i f d 3 p 1 

h = Mpf) J Wr~v 



h = 



4e 



£ _ 1 

2 



- 1 



(39) 

L x x 2 , (40) 
n. (41) 



These integrals are evaluated in Appendix [B] We obtain 
the following analytical result for the imaginary part of 
V 2 : 



ImP 2 = - 



tuAjPf) A 2 /isgn(u;) 2 / uj\ 
3T uj-Wu 2 - 4 A 2 &6C Ut/ 



|4/xw (w 2 - 4A 2 ) - T [40AV 

3w 2 (6/i + Vw 2 - 4A 2 )] sinh (^) U(|w| - 2A). 



(42) 



To leading-order the imaginary part of the polarization 
tensor is then given by 

lmn R (q)=ImP 2 (uj)(j-^ + 0(q 4 ). (43) 

The real part of the polarization tensor follows from the 
dispersion (Kramers-Kronig) relation: 



Rell^g) = / dw' 

ir J u> — u> 



(44) 



Using these quantities one can construct an effective the- 
ory of collective excitations. Their (full, interacting) 
propagators are completely determined by the spectral 
function of the collective excitations 



B(u t q) 



-2lmU R (u,q) 



[u 2 q 2 - RjeH R (u), q)] + lmU R (cj, q) 2 



The dispersion relation of the collective excitations is 
read-off as ui 2 = q 2 + ReH R (uj,q). The finite life-time 
effects are described by the width of the spectral func- 
tion, i.e., by the function Im.H R (uj,q). 

Figure Q] illustrates the dependence of the real and 
imaginary parts of the pair breaking polarization ten- 
sor on the transferred energy for fixed three-momentum 
transfer. The value of the Landau parameter is /o = 
—0.5 Q. The zero temperature gap at pf = 0.1 fm _1 
is A(0) = 1 McV and T c = A(0)/1.76. The frequen- 
cies are normalized to the zero-temperature threshold 




co/2A 



FIG. 1: Dependence of the real (solid lines) and imaginary 
(dashed lines) parts of the pair breaking polarization ten- 
sor in units of density of states on the energy transfer in 
units of threshold energy 2A(T) for fixed momentum trans- 
fer q — O.lpF, with pf = 0.1 fm _ , and two tempera- 
tures T = 0.2T C (heavy lines) and T = 0.9T C (light lines). 
The zero temperature gap is taken to be A(0) = 1. MeV, 
T c = A(0)/1.76. 



frequency 2A(0), the momentum transfer to the Fermi 
momentum. The real and imaginary parts of the polar- 
ization tensor scale as q 2 . Their behavior at negative en- 
ergies follows from their even and odd parity with respect 
to the energy transfer, i.e., ~ReH R (—uj) = KeH R (io) and 
I m II fl (— w) = — Imtl^w). Note that the imaginary parts 
are identically zero below the threshold for pair breaking 
process 2A(f). 

Figure[2]illustrates the spectral functions of pair break- 
ing density fluctuation on the energy and momentum 
transfers. The parameters are the same as in Fig.[TJ The 
form of the spectral function suggests that at low temper- 
atures the low-momentum-transfer contribution is con- 
centrated near the pair breaking threshold; for larger mo- 
mentum transfers, modes away from the energy threshold 
become important. At higher temperatures (T <T C ) and 
for any given momentum transfer, the main contribution 
to the spectral function comes from higher energy modes 
and the peak values are larger in the latter regime. 



B. Scattering response function: numerical 
evaluation 



The scattering part of the response function is kine- 
matically important for the space-like processes and van- 
ishes automatically in the limit q — > 0. Small momentum 
expansion of the previous section was found inappropri- 
ate for the scattering part of the polarization function 
and it was evaluated numerically by adapting the method 





FIG. 2: The spectral function of density fluctuations at 
T = 0.2T C (upper panel) and T = 0.9T C (lower panel). The 
parameters are as in Fig. [T] 



described in Ref . |20| • On carrying out the angular inte- 
grals we are left with a one-dimensional integral over the 
energy. As an example, we give the expression for the 
elementary loop 



lmR' GG {q)v{ PF )-\ 



2 2 
U p U p+q 



irppT 



T 



„2„,2 



e P+9 



[f(e P )-f(e p + io)]d(l- 
[f(e p )-f(e p -u)]6(l 



141) 



kol) 



(46) 



FIG. 3: Dependence of the real (solid lines) and imaginary 
(dashed lines) parts of the scattering polarization tensor in 
units of density of states on the energy transfer in units of 
threshold energy 2A(0) for fixed momentum transfer q — 
O.lpF, with pf = 0.1 fm _ , and two temperatures T = 0.2T C 
(upper panel) and T = 0.9T C (lower panel). The zero temper- 
ature gap is taken to be A(0) = 1 MeV, T c = A(0)/1.76. 



tensors Il' xx ,(q) are similar to Eq. (|46|) , but involve dif- 
ferent combinations of coherence factors. The real parts 
are obtained from the dispersion relation (|44D . 

Figure [3] shows the real and imaginary parts of the 
scattering polarization tensor at two temperatures and 
for fixed momentum transfer. Figure 2] shows the spec- 
tral function derived from the scattering polarization ten- 
sor. It is seen that the modes with energies within the 
breaking threshold u> < 2A(0) are relevant for scattering 
processes, as opposed to the pair breaking case where 
only the modes with u> > 2A(0) emerge. This justifies 
the separation of the modes into two classes. The peak 
values are larger at higher temperatures T < T c , as is the 
case for the pair breaking response. 



C. Unpaired and uncorrelated limits 

Here we wish to obtain the unpaired (T > T c ) limit of 
the polarization function [)18p. This amounts to setting 
the coherence factors in Eq. (fT9]) - (|22|) to their values in 
the normal state 



1, 



0. 



(47) 



It follows then from Eq. (|2TJ)) that Ufg = = Hgf and, 
therefore, B — 0. Setting B = in Eq. JTSfl and noting 



where 9 is the Heaviside step function, x = (p- q)/\p\\q\, that C ^ 0, since II G t G ^ 0, we obtain 



and x is the x value satisfying the equation ±u 
e p+q = an( l ^ ne P r i m e refers to the scattering part. The 
expressions for the scattering parts of the polarization 



^R.PA(q) 



A(q) 



1 - v ph A{q) ' 



(48) 
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B((o,q) 




0.1 

co/2A 



0.007 
0.006 
0.005 
0.004 
0.003 
0.002 
0.001 





0.1 

0.1 

q/p F 0.2 



0.1 

co/2A 



FIG. 4: The spectral function of density fluctuations at 
T = 0.2T C (upper panel) and T = 0.9T C (lower panel). The 
parameters are as in Fig. [T] 



where the function A(q) in the unpaired state reduces to 

d 3 P /(cp) - f(e k ) 



A(q)=Tl GG (q,A = 0)=g 



(2tt) 3 uj + e p — tk +iS 



IV. CONCLUSIONS 

In this article we have carried out several steps in the 
program aimed at understanding the polarization tensor 
of pair-correlated neutron and nuclear matter by (i) con- 
structing a low-transfered-momentum approximation to 
the pair breaking polarization tensor and (ii) by numer- 
q ically evaluating the scattering part of the polarization 
tensor. 

Our main result is that the low-momentum expan- 
sion of the pair breaking polarization tensor starts at 
quadratic order in the ratio of the momentum-transfer 
to the Fermi momentum. The expansion coefficient is 
universal, i.e., depends only the two relevant scales: the 
Fermi energy and the pairing gap. We also clarified the 
structure of the integral equations for the vertex func- 
tions when the particle-particle and particle-hole inter- 
actions do not coincide and verified explicitly that the 
unpaired and uncorrelated limits are recovered. 

Further steps will require a verification of the con- 
vergence of the series by a comparison of our analyt- 
ical results with an exact numerical evaluation of the 
response functions. The rate of series convergence can 
be checked by computing the next-to-leading-order con- 
tribution. Further refinements could include finite range 
interactions, tensor forces, and so on. The new numerical 
and analytical methods, discussed in this article, could 
be useful in the studies of the response functions of pair- 
correlated fermionic systems in general. 

The main implication of our study concerns the neu- 
trino emissivity via the vector current pair breaking 
bremsstrahlung, which are phenomenologically impor- 
tant in the physics of neutron star cooling and super- 
bursts in accreting neutron stars Our results 
suggest that the vector current contribution to the pro- 
cess 



{raj} — » {raj} + v + v, 



(51) 



(49) 



where {nn} refers to the pair correlated state, v and v to 
the neutrino and antineutrino, is suppressed to a lesser 
degree, than suggested in the recent literature 0, fl8j . 
However, further studies indicated above are needed to 
draw a final conclusion on the relative importance of this 
vector current neutrino emission process (151 



Equations (|48|) and (|49|) are the standard expressions for 
the polarization tensor of a Fermi liquid in the random 
phase approximation (RPA). The free Fermi gas result 
follows on setting v P h — in Eq. (|4"5)l : 



Ilfree(q) = 9 



d 3 P f(e p ) - f(e k ) 
(2tt) 3 u + e p — ek + i( 



(50) 



Clearly, the pair-braking contribution vanishes as T — > 
Tc . The scattering contribution at T > T c reproduces 
the unpaired and uncorrelated limits, as it should. 
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Appendix A: Expanding the functions A(q), B(q), 
and C(q) 

We write each of these functions as a product of an ap- 
propriate coherence factor and corresponding statistical 



9 



factor 



f d 3 p 

•A(q) = J -—-^a(p,q,u})L{p,q,u}), (Al) 
B(q) = Ac J -^±-L(p,q,u>), (A2) 



C(q) 

where 

a(p,q,u) 

c(p,q,u) 

L(p,q,u) 



f dp 
d 3 p 1 



(2tt) 3 2e p 



e p T e p+<J 



+ 



- 2/(e P ) 




CpCp+q 


A 2 


£p£p+q 


£p£p+q 




A 2 







(A3) 



J 



1 - /(e p ) - /(e p + g ) 
(w + iS) 2 - (e p + e p+q ) 2 



(A4) 

(A5) 
(A6) 



Next we expand these functions in powers of small pa- 
rameter y to order 0(y 2 ) 

2 A 2 2A 2 fi£ p 

a(p,q,u) = 5 — x y 

e p e p 

+ ^L(6pt p x 2 -e 2 p )y 2 , (A7) 
e P 

c(p,q,u) = 2e p + ^xy + f^ y\ (A8) 



L(p,q,u)) = L Q +L 1 xy + L 2 y 2 



(A9) 



where the coefficients of the expansion (| A9|) can be writ- 
ten using the shorthand expressions D = lu 2 — 4e 2 and 
C = e p /2T as 



L = 



tanh £ 
D + i5'' 



(A10) 



r - «„* tanhC //g P sech g , Ain 

Ll 8 ^ P JFTW + ^WT^y (An) 



4 n 2 x 2 £ p 2 sech 2 C M 2 2:2 sech 2 ( 



e p T(£ + ;<S) 2 2e p T(D + i8) 
p £ p sech 2 £ p? x 2 £ p 2 sech 2 £ 



4e p T(D + iS) 2 e p 3 T (D + iS) 
64 [i 2 x 2 i P 2 tanh ( 8 /i 2 x 2 tanh C 

(D + i5f (D + iSf 

4 /I £ p tanh £ 4 p? x 2 £ p 2 tanh ( 
(D + iS) 2 e 2 (D + i8f 

,2Jt2 



jj, x £ p tanh £ sech £ 
2e p 2 T 2 (D + z<5) ' 



(A12) 



where sech 2 (x) = 1 — tanh(x) 2 . On substituting 
Eqs. (TAT)) and flAjj) in Eq. (JATJ) we obtain 

~2A 2 , 2A 2 



2A 2 M £ P 



-La 



T 2 2 

Li x y 



-L 2 y A 

t P 

A 2 m^p, 



5 (6^£ P x -e p )L y 

2 i r>f„.4 



A>(<z) + -4 2 (<zh/ 2 + 0(2/ 4 ) 



(A13) 



The terms that are odd in x drop out on integration in 
symmetrical limits. Similarly, 



C(q) 



d 3 p 
(2Vp 

2e p L 2 y 2 



2e p L a + 2 -^L x x 2 y 2 + ^ L y 2 



d 3 p 1 



(2tt) 3 2e p 
We further use the relation 



[1 - 2f(e p )} . (A14) 



d 3 p 
(2^)3 



2e p Lo 



d 3 p 1 

JW 3 ^ 



[1 - 2/(e p )] 
d 3 p 1 



2 ./ ( 27 r)3 £p L0 (A15) 



to write 



C(q) = - 



d 3 p 1 
2 ./ (2^)3"^ 



d 3 p 



2^Cp r 2 2 

-Li x y 



^ L y 2 + 2e p L 2 y 2 

673 



(2tt) 3 
--C (q)+C 2 (q)y 2 



(A16) 



Finally, 

B(q) = Auj f J^-±(Lo + L 2 y 2 ) = B (q)+B 2 (q)y 2 , 
J (Zir) a e p 

(A17) 

and Ao/Co = -4A 2 /cj 2 , and B /C = -2A/w. 

Appendix B: Evaluating phase-space integrals 

Here we evaluate the integrals (|59^ -(f4*T j) . which are 
given as 



Io = v{pf) V 



d 3 p £ p 



(2^) 3 e 3 



(Bl) 



4e; 



(2tt) 3 e F 



4e 2 
p _ 1 

w 2 



Lix 2 , (B2) 
£2(1). (B3) 
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After carrying out the angular integrals we obtain For positive lo the integral (jBll) gives [2l| 



Io = vipp) fi J 

2 i fdpp^p 



K 1 2/xg p 



4e 2 



J 2 = 2z/(j?i^ 



-l 



dpp 2 1 
2tt 2 Z 



4e 2 
c p _ i 

LO 2 



(B4) 
£i, (B5) 

(La(i)), (B6) 



where (L2) is the angle integrated loop L2, which is ob- 
tain from (|A12[) via the substitution x 2 — 1/3. Using the 
transformation dpp 2 = m*ppd{p 2 /2m*) = m*ppd£ l p, and 
the relation e p de p — we obtain 



Io = M 



dCp^r 



4e 2 
p _ 1 

w 2 



r 2 ^ 1 



4e 2 
c p _ 1 

2 



I2 = 2 



?p 



.2 

2 - 1 
,2 



io, (B7) 
Li, (B8) 
(L 2 (x)). (B9) 



To make further progress we need to separate the real and 
imaginary parts of the integrals. We shall first compute 
the imaginary parts. They are extracted with the help 
of the identity 

1 P (—'\') n 

77^ 777— —7 = "7^~~tt - iir^-r-SW (£>), (BIO) 

(D + i5) n+1 D n+1 n\ 

where P denotes the principal value and S^ n \D) is the 
ri-th derivative of the delta function. 



ImJo 



-7T/1 



dc 



P £ 2 

p 



tanh(^)^ 2 -4 e 2 ) 



o , Vw 2 - 4A 2 sgn(w) , / lo 



{lo - 2A). 
(Bll) 



Consider the integral (IB2[) . First, note that the term oc 
D^ 1 vanishes, since the delta function enforces lo 2 = 4e 2 . 
After dropping this term we are left with the integral 



Im/i 



(B12) 



This and similar integrals, which contain derivatives of 
the delta function, are computed via the formula 



f(x)8 n (x-a)dx = (-) n f( n \a). 
The result of integration is 



(B13) 



ImJx 



8V V^-4A 2 sgn H tanh ^ 9{u _ 2Ay 

3 lo 5 \4TJ y ' 

(B14) 

In the integral (|B3I) we again omit terms « D 1 , since 
their prefactors are zero after integration. After inserting 
x 2 = 1/3 in the remainder we obtain 



Im/2 = — I de 

A 



4e 2 
v _ 1 

lo 2 



2 J el - A 2 



e v T 



-sech 2 ^ 
2T 



e 2 - A 2 



+3/^tanh(^) 



» 2 Je 2 p -A 2 



tanh 



u 2 tanh 

P V2T 



64tt^ 2 



2 /-oo 



A 



4et 



e 2 -A 2 tanh(|l)^V-4 e 2 ). 



(B15) 



Applying Eq. (|B13|) we obtain 
■Kji sgn(w) 



Mo 



sech' 



(JjL^ |r [24AV + lo 2 (2/i + 3^ 2 -4A 2 )] sinh (^) + Afiu {lo 2 - 4A 2 )| 6»(o; - 2A) . 



3Tlo 5 Vlo 2 - 4A 2 

Finally, adding the three integrals (jBllI) . (|B14|) . and (|B16|) we obtain Eq. (|4"2")l of the main text. 



(B16) 
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